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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  at  the  request  of  the  AEDC/DOT. 
The  AEDC  Project  Manager  was  Mr.  A.  F.  Money.  The  results  of  the  research  were 
obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup  Corporation  Company),  operating 
contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO  Project 
Number  P32E-39D.  The  manuscript  was  submitted  for  publication  on  October  8,  1979. 

This  report  is  the  first  in  a  series  of  four  volumes  entitled  “Store  Separation  Testing 
Techniques  at  the  Arnold  Engineering  Development  Center."  Subtitles  of  these  volumes  are 
as  follows: 


Volume  1 

An  Overview 

Volume  11 

Description  of  Captive  Trajectory  Store  Separation  Testing 
in  the  Aerodynamic  Wind  Tunnel  (4T) 

Volume  111 

Description  and  Validation  of  Captive  Trajectory  Store 
Separation  Testing  in  the  von  Karman  Facility 

Volume  IV 

Description  of  Dynamic  Drop  Store  Separation  Testing 
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1.0  INTRODUCTION 

As  the  operational  flight  speeds  of  military  aircraft  have  increased,  the  problems 
associated  with  aircraft/store  compatibility  have  increased  significantly.  Generally,  this  can 
be  attributed  to  a  large  increase  in  the  ratio  of  store  aerodynamic  to  inertial  forces  caused  by 
the  increased  flight  speeds  and  to  independent  rather  than  integrated  design  of  aircraft  and 
stores.  It  is  obvious  that  compatibility  for  the  large  number  of  aircraft/store  combinations 
cannot  be  determined  from  full-scale  flight  tests  alone:  not  only  would  the  costs  be  pro¬ 
hibitive;  the  unnecessary  risks  to  the  pilot  would  also  be  unacceptable.  Therefore,  certifica¬ 
tion  of  aircraft/store  load  combinations  for  carriage  and  release  must  be  obtained  by  other 
means  and  confirmed,  as  required,  by  flight  tests. 

Currently,  many  facets  of  aircraft/store  compatibility  are  evaluated  and  verified  from 
wind  tunnel  test  data  obtained  on  reduced  scale  models.  For  instance,  aircraft  stability,  car¬ 
riage  loads,  captive  separation  trajectory,  aerodynamic  grid  survey,  flow-field  survey,  store 
heat  transfer,  and  dynamic  drop  separation  trajectory  tests  are  routinely  conducted  in  the 
AEDC  wind  tunnels.  However,  this  report  considers  only  separation  trajectory  evaluation. 

A  separation  trajectory  can  be  defined  as  a  time  history  of  the  position  and  attitude  of  a 
released  store  with  respect  to  the  aircraft.  Separation  trajectories  generally  fall  into  three 
categories.  First,  the  released  store  may  undergo  a  clean  separation  from  the  aircraft  in  a 
manner  that  allows  a  high  probability  of  striking  the  desired  target.  Second,  the  released 
store  may  separate  cleanly,  but  an  unfavorable  aircraft  flow  field  could  induce  large  pitch 
and/or  yaw  excursions  such  that  on-board  guidance  equipment  would  fail  to  lock  on  target 
(guided  weapon)  or  could  so  alter  the  initial  ballistics  of  an  unguided  weapon  as  to  cause 
impact  well  outside  the  acceptable  targeting  limits.  Third,  the  aircraft  flow  field  could 
conceivably  cause  the  released  store  to  “fly”  back  into  the  aircraft  with  potentially 
disastrous  results.  Obviously,  no  store  would  be  certified  for  release  at  conditions  where  any 
real  possibility  of  an  aircraft/store  collision  existed.  It  would  also  be  undesirable  to  release 
munitions  under  conditions  that  precluded  any  reasonable  chance  of  a  successful  strike, 
even  though  the  separations  might  be  clean.  Therefore,  it  is  highly  desirable,  if  not 
imperative,  to  define  release  conditions  where  the  risk  is  lowest  and  the  success  probability 
greatest  before  flight  testing  is  initiated. 

Three  basic  approaches  are  available  for  generating  separation  trajectories:  dynamic 
drop  wind  tunnel  tests,  wind  tunnel  computer  simulations,  and  analytic  computer 
simulations.  The  dynamic  drop  approach  is  strictly  experimental  while  the  computer 


7 


AEDC-TR-79-1 


simulations  may  use  experimental,  semiempirical,  or  theoretical  techniques  (see  block 
diagram,  Fig.  1).  The  purpose  of  this  report  is  to  describe  several  of  the  options  by  which 
store  separation  characteristics  may  be  generated  and  discuss  the  relative  merits  and 
shortcomings  of  each  such  that  guidelines  for  their  general  use  may  be  implied.  This  report  is 
the  first  in  a  series  of  four  volumes  which  describe  store  separation  capabilities  at  the  AEDC. 
Volumes  II  and  III  describe  captive  trajectory  store  separation  testing  in  the  AEDC 
Aerodynamic  Wind  Tunnel  (4T)  and  Supersonic  Wind  Tunnel  (A),  respectively.  Dynamic 
drop  testing  capabilities  for  all  AEDC  wind  tunnels  are  given  in  Volume  IV. 

2.0  DYNAMIC  DROP  TRAJECTORY  TECHNIQUES 


2.1  GENERAL 

As  previously  noted,  separation  trajectories  by  the  dynamic  drop  (or  free-drop)  method 
must  be  effected  experimentally.  However,  because  of  the  nonsteady-state  nature  of  the 
tests,  it  is  necessary  to  simulate  not  only  the  applied  forces,  but  also  the  inertial  response  of 
the  store  to  these  forces  (Ref.  1).  The  aerodynamic  forces  can  be  simulated  on  a 
geometrically  scaled  model  if  the  compressibility  and  viscous  airstream  effects  are 
reproduced.  The  principles  of  dynamic  similitude  for  reduced  scale  models  are  not  new  and 
have  been  well  documented  (e.g.,  Refs.  1  through  3).  However,  the  method  shown  by 
Marshall  (Ref.  1)  will  be  used  to  establish  the  dynamic  relationships  for  this  report. 

2.2  THREE-DOE  EQUATIONS  OF  MOTION 

By  using  simplified  (three  degrees  of  freedom)  equations  of  motion  and  an  earth- 
oriented  coordinate  system  (see  Fig.  2),  the  forces  and  moments  on  the  store  are  described 
as: 


mX  =  — 

|~Ca  t  cose 
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— 
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II 
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1  ' 

V 
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(3) 
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r  El  aE  1 

If  the  angular  motion  is  small  and  if  flow-field  spatial  variations  in  streamline  inclination 
and  curvature  can  be  represented  by  incremental  values  of  angle  of  attack  and  pitching- 
moment  coefficient,  respectively,  then  the  aerodynamic  coefficients  may  be  described  as: 

C»-CN.(*^+A*0  (4) 
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c  =  c  (0+A+Aa)+AC 
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After  the  coefficient  relationships  are  substituted,  the  motion  equations  become: 
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2.3  DYNAMIC  SCALING  PRINCIPLES 


(5) 


(6) 


(7) 


(8) 


(9) 


From  an  examination  of  the  previous  three  equations,  one  can  deduce  the  dynamic 
simulation  relationships  which  will  produce  comparable  motion  on  model  and  full-scale 
stores.  From  the  reejuirernent  of  linear  geometric  scaling  of  stores  and  aircraft  models  (where 
X  is  the  scale  factor),  the  relationships  are  as  follows: 

1 


A'  = 

AX2 

e'  =  e 

*1- 

EjA 

X'  =  X  A 

eg  eg 
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XA 

XEl  =  XElX 
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where  primed  quantities  represent  model  quantities. 


If  compressibility  and  viscous  effects  in  the  airstream  are  matched  on  the  geometrically 
scaled  models,  then: 


CN 

=  CN 

C  A,t 

=  CA,t 

a 

a 

C' 

=  C 

C' 

=  c 

m  m 

q  q 


(id 


From  the  requirement  that  full-  and  model-scale  values  of  (X/g),  (Z/g)  and  8  be  equal,  the 
full-  and  model-scale  values  of  the  parameters  (Z/Voo),  (qooA/mg),  (Fei/mg),  (f]0/2Voo), 
(qooA  fi/I  yy)j  and  (Fei  XEi/Iyy)  must  also  be  equal.  The  relationships  of  time,  velocity  (and 
Mach  number),  mass,  inertia,  and  ejector  force,  respectively,  are  then  determined  to  be: 


t'  =  WA(g/gO 

(12) 

V'oo  Voc  VA( g'/g) 

(13) 

V  A(gVg)  (T^/T^) 

(14) 

=  g/g0A2 

(15) 

lYY  =  1yyWpJ  (V^/VJ2(g/g')A4 

(16) 

fei  =  fei  (p:/W(to2a2 

(17) 

Noting  that  g'  =  g  and  substituting  the  velocity  relationship  into  the  other  equations  gives 


t'  = 

t  v/  A. 

(18) 

V'  = 

CO 

V  \TT 

OO  7 

(19) 

M'  = 

oo 

(20) 
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(21) 

lYY  = 
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Equations  (18)  through  (23)  are  generally  referred  to  as  the  “Froude  scaling”  relationships 
because  the  velocity  scaling  is  equivalent  to  the  hydrodynamic  Froude  number. 

As  noted,  the  only  stipulations  for  the  Froude  scaling  principles  are  linear  geometric 
scaling  and  aerodynamic  coefficient  equality.  These,  in  turn,  lead  to  the  reduced  velocity 
requirement  (assuming  g'  =  g)  on  the  scaled  models.  As  long  as  full-scale  flight  conditions 
are  low  subsonic,  these  relationships  are  compatible  and  valid.  However,  with  the  onset  of 
significant  airstream  compressibility  effects,  model-scale  aerodynamic  data  obtained  at 
reduced  velocity  may  no  longer  be  equal  to  the  full-scale  values.  Therefore,  the  useful  test 
range  of  the  Froude  relationships  is  restricted  to  low  subsonic  flight  conditions. 

To  insure  aerodynamics  matching,  compromise  simulation  methods  have  been 
developed  which  require  Mach  number  equality  at  the  expense  of  another  parameter.  These 
Mach  scaling  techniques  are  commonly  known  as  the  “heavy”  and  “light”  scaling 
principles.  For  heavy  scaling,  which  is  probably  the  most  widely  used  scaling  technique,  the 
quantities  {Z/V  m)  and  (f)0/2Voo)  are  removed  from  the  simulation.  The  relationships  for 
heavy  scaling  which  are  changed  from  Froude  scaling  are  then  (assuming  g'  =  g): 


M'  =  M 

DO  OQ 

(24) 

V'  =  V  \l  ('!”/'!'  ) 

OO  DO  V  oo'  do' 

(25) 

m'  = 

(26) 

YY  =  ^Y^/O*4 

(27) 

fei  =  FEi(<w/qJ;'-2 

(28) 

The  connotation  “heavy  scaling”  may  be  deduced  from  the  mass  relationship.  All  else  being 
equal,  the  mass  of  the  heavy-scaled  models  will  be  larger  than  that  for  Froude-scaled  models 
by  a  factor  of  (1/A)  because  of  the  increased  velocity  requirement.  Because  of  the  relaxed 
velocity  ratio  simulation,  the  heavy  scaling  principles  do  not  properly  account  for  induced 
angle  of  attack  or  aerodynamic  damping  effects.  The  amplitude  of  the  angular  motions  will 
be  too  large  because  of  reduced  damping,  and  linear  position  will  be  affected  by  the 
dependence  of  the  aerodynamic  coefficients  on  angle  of  attack.  However,  previous  work 
(Ref.  2)  has  indicated  that  the  angle-of-attack  error  is  never  greater  than  10  percent  and 
that  the  damping  error  is  of  small  consequence  in  studying  body  motion  immediately  after 
release,  unless  the  rate  of  decay  of  oscillation  is  important. 
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When  the  angular  motion  response  is  deemed  important,  the  light  scaling  principles  may 
be  used.  Light  scaling  retains  the  velocity  ratio  simulations  by  assuming  that  the 
gravitational  constant  of  the  model-scale  experiment  may  be  arbitrarily  increased,  such  that: 

s'  ‘  &)'/'  <»> 
The  scaling  relationships  for  light  scaling  which  are  changed  from  Froude  scaling  are  then: 

M'  =  M 

oo  oo 

=  vJt;/tm 

FE1  = 

t'  =  {viv./vj] 


(30) 

(31) 

(32) 

(33) 


The  connotation  “light  scaling”  is  derived  from  the  equality  of  the  mass  ratio  with  Froude 
scaling.  Unless  the  increased  gravitational  requirement  can  be  effected  with  a  magnetic  field 
(Ref.  4)  or  perhaps  by  accelerating  the  aircraft  away  from  the  store  after  separation  occurs, 
the  model  will,  in  reality,  be  subjected  to  the  gravitational  acceleration  “g”  rather  than 
“g'”.  The  vertical  acceleration  will  then  be  deficient  by  the  amount 


Ag'  = 


■(V'/V  )2 

x  no'  no 


(34) 


Consequently,  vertical  displacement  will  be  too  small.  Up  to  the  end  of  the  ejector  stroke, 
the  correct  motion  can  be  produced  with  an  increased  ejector  force  such  that: 


afei  =  m'Ag" 


(35) 


XE1  =  XE1A  FEi^FE1  +  AFEl|] 


(36) 


Even  so,  with  the  gravitational  acceleration  deficient  throughout  the  remainder  of  the 
trajectory,  the  risk  of  premature  collision  is  increased.  To  avoid  this  possibility,  the  ejector 
force  could  be  further  augmented  such  that  a  velocity  simulation  match  would  be  effected  at 
a  point  in  the  flow  field  other  than  the  end  of  the  ejector  stroke  (typically  one-half  store 
length  away  from  carriage).  However,  the  velocity  ratio  simulation  would  be  lost  on  either 
side  of  the  match  point. 
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Marshall  (Ref.  1)  has  shown  the  feasibility  of  the  augmented  ejector  force  concept  for 
light  scaling  and  also  has  developed  a  correction  factor  which  may  be  applied  to  the  model- 
scale  trajectory  to  account  for  insufficient  gravity  and  the  increased  ejector  force.  The 
correction  is  an  empirical  relation  which  adjusts  the  model-scale  vertical  displacement  curve 
and  is  defined  as: 


/AF'  \ 

AZ'=  V2(Ag')  (to2  -  Va(— ^  )  (to2  for(t'<t£) 

\  m  / 


(37) 


az'  =  wxt')2-  y) 

\  m  V 


(tf)(2t'-tp)  for  (t'>  tA) 


(38) 


As  may  be  seen,  use  of  the  correction  requires  that  the  time  at  the  end  of  the  ejector  stroke 
(t^)  be  known.  If,  however,  the  value  of  (tp)  is  unknown,  it  can  be  estimated  if  one 
assumes  no  aerodynamic  forces  during  ejector  action.  The  resulting  error  in  the  corrected 
displacements  typically  will  be  on  the  order  of  2  percent  (Ref.  1). 

2.4  TRAJECTORY  COMPARISONS 


To  assess  the  relative  merits  and  disadvantages  of  the  dynamic  scaling  techniques,  direct 
comparison  of  trajectories  by  each  of  the  techniques  with  full-scale  flight  test  trajectory 
results  would  be  desirable.  However,  in  the  absence  of  flight  test  data,  analytical  trajectories 
were  generated  using  the  3-DOF  equations  of  motion  of  Section  2.2  (Ref.  1).  (For  reference, 
the  Froude,  heavy,  and  light  scaling  principles  are  summarized  in  Fig.  3).  For  full-scale  flight 
conditions  of  =  0.95  and  altitude  of  6,000  ft,  trajectories  were  calculated  on  a  2,000-lb 
store  (large,  slender,  finned  bomb)  using  aerodynamic  data  as  typically  expected  for  the 
store  and  aircraft  flow  field  (see  Figs.  4  and  5,  respectively).  Using  the  same  motion 
equations  and  aerodynamic  data  (except  as  noted),  trajectories  were  then  calculated  for 
wind  tunnel  models  which  were  designed  according  to  each  set  of  scaling  laws  and  for  the 
following  conditions: 


0.05 

=  50  psf  (Froude) 

570°R 

=  200  psf  (Heavy) 

q^  =  600  psf  (Light) 

These  trajectory  results  are  presented  in  Figs.  6  through  8.  When  the  transonic  aerodynamic 
coefficient  data  of  the  full-scale  store  were  used  with  the  Froude-scaled  model  (Fig.  6a),  the 
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model  motion  reproduced  that  of  the  full-scale  store  exactly.  However,  using  transonic 
aerodynamic  characteristics  at  the  reduced  velocity  required  for  Froude  scaling  effectively 
ignores  changes  in  aerodynamics  between  compressible  and  incompressible  flow.  When 
typical  subsonic  aerodynamic  data  are  used  with  the  Froude-scaled  models  (See  Figs.  4  and 
5),  differences  in  full-  and  model-scale  trajectories  are  large  (Fig.  6b). 

The  heavy  scaling  trajectory  data  indicate  an  underdamped  pitch  oscillation  for  the 
model-scale  trajectory  because  of  incorrect  velocity  scaling  (Fig.  7).  However,  the  agreement 
obtained  between  full-  and  model-scale  trajectory  data  is  quite  reasonable  for  the  angular 
motion  and  quite  good  for  the  linear  motions.  The  conclusions  drawn  from  the  heavy-scaled 
analytical  data  are  confirmed  by  wind  tunnel  and  flight  test  results  from  the  F-15  aircraft 
(see  Fig.  9  and  Ref.  5). 

Calculated  trajectories  for  light-scaled  models  were  determined  for  three  conditions:  (1) 
no  ejector  augmentation  (Fig.  8a),  (2)  ejector  force  increased  by  the  amount  necessary  to 
offset  gravitational  deficiency  during  the  ejector  stroke  (see  Section  2.3  and  Fig.  8b),  and  (3) 
ejector  force  augmented  by  three  times  the  amount  used  in  the  preceding  case  which 
provided  the  velocity  simulation  match  near  one-half  body  length  from  carriage  (Fig.  8c). 
For  all  three  cases,  the  angular  motion  of  the  full-scale  trajectory  was  well  matched  by  the 
model  data,  but,  as  expected,  vertical  position  was  considerably  in  error.  However,  when  the 
vertical  position  correction  factor  (Ref.  1)  was  applied  to  the  model  data,  agreement  was 
quite  reasonable.  It  should  be  noted  that  for  all  three  scaling  principles,  similar  trends  were 
exhibited  in  trajectory  data  calculated  for  an  unstable  store  of  the  750-lb  class  (Ref.  1)  such 
that  the  stable  store  results  are  not  considered  unique. 

2.5  APPLICATIONS 

The  choice  among  the  available  techniques  for  the  design  of  dynamic  drop  models  is  not 
an  easy  one.  From  a  purely  analytical  standpoint,  the  selection  of  Froude  scaling  would  be 
desirable  since  full-scale  motion  is  reproduced  exactly  within  the  limits  of  the  simulation. 
However,  as  previously  mentioned,  to  choose  Froude  scaling  with  confidence  would  require 
almost  certain  knowledge  that  the  store  aerodynamic  coefficients  were  invariant  with  Mach 
number.  Unfortunately,  experience  has  shown  that  aerodynamic  data  obtained  around 
typical  carriage  locations  at  transonic  and  supersonic  speeds  will  almost  certainly  be  Mach 
number  dependent;  this  condition  effectively  eliminates  the  choice  of  Froude  scaling  for 
many  applications. 

The  trajectory  data  of  Figs.  7  and  9  seem  to  substantiate  the  usual  choice  of  heavy  scaling 
as  the  alternative  to  Froude  scaling.  In  all  but  the  most  demanding  analyses, 
expected  differences  between  model-  and  full-scale  trajectory  data  would  surely  be  within 
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acceptable  limits.  However,  there  are  some  practical  problems  inherent  in  this  technique. 
The  large  model  mass  requirement  quite  often  necessitates  the  use  of  high  density  materials 
such  as  tungsten,  gold,  iridium,  or  platinum  in  model  fabrication.  Disadvantages  of  these 
materials  generally  include  one  or  more  of  the  following:  they  arc  1)  not  readily  available,  2) 
difficult  to  machine,  or  3)  very  expensive.  Reduction  of  the  dynamic  pressure  for  the  wind 
tunnel  experiment  could  alleviate  the  high  mass  requirement  somewhat  but  only  at  the 
further  expense  of  the  Reynolds  number  simulation. 

A  viable  scaling  alternative  is  the  light  scaling  technique  since  it  has  been  demonstrated 
(Ref.  1)  that  position  error  from  insufficient  vertical  acceleration  can  be  corrected. 
However,  to  establish  the  proper  experimental  parameters,  a  general  knowledge  of  store 
behavior  must  be  obtained  in  advance.  Unfortunately,  one  of  the  primary  reasons  for  free- 
drop  testing  is  that  a  nonuniform  and  largely  undefined  aerodynamic  flow  field  surrounds 
an  aircraft  at  transonic  and  supersonic  speeds.  An  additional  disadvantage  of  this  technique 
is  that  implementation  of  a  large  augmented  ejector  force  on  small-scale  models  may 
generate  a  completely  new  set  of  design  problems  which  are  Comparable  in  complexity  to 
those  of  heavy  scaling. 

The  previous  arguments  have  been  confined  to  the  relative  merits  and  shortcomings  of 
each  of  the  dynamic  drop  techniques  when  compared  to  the  others.  However,  the 
advantages  and  disadvantages  of  the  dynamic  drop  experiment  with  respect  to  other 
trajectory  acquisition  techniques  must  also  be  considered  in  any  practical  applications.  Some 
of  the  characteristics  of  the  dynamic  drop  technique  which  may  be  considered  advantageous 
are  as  follows: 

1.  Test  results  are  valid  within  the  limits  of  the  scaling  relationships. 

2.  Since  no  modifications  are  required  on  the  reduced-scale  stores  to  accommodate 
balance/sling  combinations,  no  errors  in  the  trajectory  data  will  be  introduced 
from  anomalies  in  aerodynamic  data  caused  by  either  configurational 
compromise  or  sting  interference. 

3.  For  weapons  bay  releases,  trajectories  may  be  initiated  from  the 
carriage  position. 

4.  There  are  no  hardware  restraints  to  obstruct  the  motion  of  tumbling  stores. 

5.  Several  types  of  staged  release  separations  can  be  simulated  (for  instance, 
pivoting  fuel  tank). 
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6.  The  dynamic  drop  technique  is  particularly  suited  for  simulation  of  complicated 
loading  separations  such  as  multiple  racks  with  a  full  or  partial  complement 
of  stores. 

7.  Multiple  or  ripple  releases  can  be  accomplished. 

Some  of  the  characteristics  of  the  dynamic  drop  technique  which  might  be  considered 
disadvantageous  are  as  follows: 

1.  The  use  of  a  wind  tunnel  facility  is  required,  and  facility  productivity  must  be 
considered  in  resource  management.  For  the  AEDC  Aerodynamic  Wind  Tunnel 
(4T)  and  Propulsion  Wind  Tunnel  ( 1ST),  approximately  1-1  /2  to  2  drops  per  test 
hour  can  be  accomplished  for  a  single  drop  mode.  Because  of  the  model 
injection  system  of  the  AEDC  Supersonic  Wind  Tunnel  (A),  it  is  estimated  that  3 
to  4  drops  per  test  hour  could  be  accomplished.  (When  multiple  drops  can  be 
made  during  one  run,  productivity  is  correspondingly  increased  for  all 

'  l 

AEDC  tunnels.) 

2.  A  separate  model  is  normally  required  for  each  trajectory.  Although  some 
models  may  be  captured  in  a  reusable  state,  model-catching  screens  are  generally 
used  only  to  salvage  models  for  refurbishing  and  use  in  subsequent 
test  programs. 

3.  Some  store  detail  may  be  lost  at  reduced  scale  (for  instance,  lugs,  conduits, 
radomes,  etc.). 

4.  Models  must  be  designed  for  the  capabilities  of  the  test  facility.  (This 
requirement  is  less  critical  for  a  variable  density  wind  tunnel  since  a  range  of 
model  weight  would  be  valid  for  a  given  set  of  flight  conditions)  (see  Fig.  3). 

5.  In  the  design  and  fabrication  of  wind  tunnel  models,  it  is  often  difficult  to  obtain 
the  correct  model  inertia,  weight,  and  center  of  gravity  location  simultaneously. 

6.  Only  straight  and  level  flight  can  be  simulated.  However,  in  much  the  same 
manner  that  vertical  position  corrections  for  insufficient  gravity  can  be  made  for 
light  scaling  (see  Section  2.3),  corrections  could  be  applied  to  the  axial  and 
vertical  position  data  for  first-order  effects  of  climbing  or  diving  flight. 
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7.  Thrusting  or  active  guidance  system  releases  cannot  be  accomplished.  Only 
jettison  trajectories  may  be  simulated. 

8.  Because  of  processing  and  reading  requirements  of  the  data  film,  access  to  final 
trajectory  data  is  comparatively  slow.  Therefore,  any  decisions  required  to 
define  test  program  direction  must  be  based  on  visual  interpretation  of  fast- 
processed  film  data. 

To  summarize  briefly,  the  overall  objectives  of  any  separation  analysis  will  dictate 
whether  the  acquisition  of  dynamic  drop  trajectory  data  is  desirable.  When  Tree-drop 
trajectories  are  deemed  necessary,  selection  of  the  proper  scaling  technique  generally 
requires  compromises.  The  choice  which  will  be  most  advantageous  to  the  particular 
experiment  can  be  made  only  by  familiarity  with  all  the  available  options.  A  complete 
description  of  how  the  dynamic  drop  experimental  data  are  obtained  and  evaluated  at  the 
AEDC  is  given  in  Volume  IV  of  this  series. 

3.0  COMPUTER  SIMULATION  TRAJECTORY  TECHNIQUES 
3.1  GENERAL 

Separation  trajectories  can  be  calculated  on  high-speed,  digital  computers  if  the  full-scale 
forces  and  moments  at  any  point  along  the  trajectory  path  and  the  equations  of  motion  of 
the  released  stores  are  known.  A  simplified  block  diagram  showing  essential  elements  of  a 
trajectory  generation  computer  program  is  presented  in  Fig.  10.  The  composition  within  the 
blocks  may  vary  considerably  among  different  trajectory  generation  programs  depending  on 
the  complexity  of  1  he  specific  program,  but  all  trajectory  programs  must  perform  in  some 
manner  the  functions  indicated. 

Although  the  equations  of  motion  are  universal,  the  boundary  conditions  imposed 
dictate  which  terms  must  be  retained  for  the  simulation  (free  or  restrained  motion,  number 
of  degrees  of  freedom  desired,  etc.).  In  addition  to  the  static  aerodynamic  forces  and 
moments,  the  full-scale  forces  and  moments  in  a  comprehensive  program  will  include  the 
effects  of  aerodynamic  damping,  gravity,  ejector  forces  and  moments,  thrust,  jet  damping, 
reaction  forces  and  moments  (to  simulate  staged  release),  and  control  inputs  (to  simulate 
active  guidance  systems,  flap  deflections,  drogue  chute  deployment,  etc.).  The  accelerations 
and  velocities  of  the  store  can  be  numerically  integrated  using  any  of  several  available 
methods  (Ref.  6).  However,  the  greatest  variations  in  computer  simulation  trajectory 
techniques  probably  exist  in  the  method  of  defining  aerodynamic  coefficient  data.  Basically, 
two  opiions  are  available  (see  Block  diagram.  Fig.  I).  At  each  point  along  the  trajectory 
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path,  the  aerodynamic  coefficient  data  can  be  measured  in  the  wind  tunnel  using  reduced- 
scale  models  and  a  captive  trajectory  support  (CTS)  system  or  predicted  using  a  flow-field 
math  model.  Each  method  will  be  described  more  fully  in  subsequent  sections  of  this  report. 

3.2  AEDC  TRAJECTORY  GENERATION  PROGRAM 

At  the  AEDC,  a  basic  trajectory  generation  package  has  been  developed  which  can  be 
used  in  conjunction  with  either  the  captive  trajectory  support  systems  or  flow-field  math 
models  for  separation  trajectory  generation.  The  general  composition  of  this  package  is 
outlined  as  follows.  The  equations  of  motion  allow  for  six-degrees-of-freedom  movement  of 
the  released  store  and  are  given  in  Appendixes  A,  B,  and  C  (see  Fig.  1 1  for  definition  of  the 
body-axis  system).  In  addition  to  free  motion  releases,  the  equations  include  provisions  for 
staged  separation  (e.g.,  pivoting  fuel  tanks,  missile  rail  launches)  and  for  accelerated  flight 
(pullup/pushovcr  aircraft  maneuver  only,  see  Fig.  12).  Aircraft  dive  or  bank  maneuvers  are 
simulated  by  adjusting  the  weight  components  of  the  store  at  release.  Full-scale  force  and 
moment  equations  for  a  free-falling  store  are  presented  in  Appendix  D.  As  may  be  seen, 
terms  are  included  to.accouni  for  each  of  the  items  mentioned  in  Section  3.1.  Of  the  force 
and  moment  parameters  included,  only  the  aerodynamic  damping  derivatives  are  normally 
assumed  to  be  constant  throughout  a  trajectory.  Integration  of  the  accelerations  and 
velocities  is  accomplished  using  the  Adams-Moulton  (predict  or -corrector)  algorithm  with  a 
Runge-Kutta  algorithm  to  start  the  process  (Ref.  6).  To  permit  the  trajectory  calculations  in 
the  corrector  loop  of  the  integration  process  to  be  made  independently  of  data  acquisition, 
aerodynamic  coefficient  data  (CY,  CN,  etc.)  used  in  the  force  and  moment  calculations  are 
obtained  by  extrapolation.  The  extrapolation  is  accomplished  using  a  second-order 
polynomial  to  fit  the  time  history  of  the  last  three  calculated  values  of  the  coefficients.  A 
complete  description  of  the  trajectory  program  is  contained  in  Volume  II. of  this  series. 

3.3  WIND  TUNNEL  APPLICATIONS 

Much  of  the  initial  work  in  store  separation  testing  was  performed  by  the  Naval  Ship 
Research  and  Development  Center  (NSRDC)  (Refs.  7  and  8).  Some  of  the  other  wind  tunnel 
facilities  which  conduct  store  separation  tests  include  those  of  the  LTV  Aerospace 
Corporation  (Ref.  9),  the  Office  National  d’Etudes  et  de  Recherches  Aerospatiales 
(ONERA,  Ref.  10),  the  Weapons  Systems  Research  Laboratory  (WSRL,  Ref.  11),  the 
Cornell  Laboratory  (Cal span,  Ref.  12),  and  the  Aircraft  Research  Association,  LTD.  (ARA, 
Ref.  13). 

For  store  separation  testing  using  computer  simulations,  (wo  separate  and  independent 
support  systems  are  used  to  support  the  models.  The  aircraft  mode!  is  attached  to  the  tunnel 
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main  support,  and  the  store  model  is  attached  to  the  captive  trajectory  support  (CTS) 
mechanism,  which  provides  store  movement  (generally  six  degrees  of  freedom)  independent 
of  the  aircraft  model.  An  isometric  drawing  of  a  typical  store  installation  in  Tunnel  4T  is 
shown  in  Fig.  13.  The  CTS  mechanism  for  Tunnel  4T  has  been  operational  since  1968,  and 
the  initial  description  of  system  operation  (Ref.  14)  is  superseded  by  Volume  II  of  this  scries. 
The  CTS  mechanism  for  AEDC  Tunnels  A,  B,  and  C  (Fig.  14)  has  been  operational  since 
1973,  and  a  description  of  system  operation  is  contained  in  Volume  111  of  this  series. 

Wind  tunnel  tests  which  utilize  either  AEDC  CTS  system  can  be  conducted  in  either  the 
trajectory  generation  or  grid  mode.  In  the  trajectory  mode,  the  trajectory  generation 
computer  program  (Section  3.2)  is  coupled  with  the  CTS  system  through  the  facility  digital 
computer  in  a  closed-loop  operation  (Figs.  10  and  13).  After  the  wind  tunnel  test  conditions 
are  established,  aircraft  model  angle  of  attack  is  set,  and  the  store  model  is  located  at  the 
carriage  position,  the  trajectory  generation  process  is  initialed.  Each  point  along  the 
trajectory  is  calculated  using  measured  aerodynamic  coefficient  data  and  the  applicable 
simulation  parameters  as  inputs  to  the  trajectory  program.  As  each  succeeding  point  is 
determined,  the  store  model  is  automatically  repositioned  through  commands  to  the  CTS, 
and  new  aerodynamic  data  inputs  are  obtained.  This  process  is  repeated  until  trajectory 
completion.  As  may  be  seen,  the  only  essential  inputs  provided  by  the  wind  tunnel  for  the 
simulation  are  the  aerodynamic  coefficient  data.  Typical  trajectories  as  obtained  in  Tunnel 
4T  are  presented  in  Fig.  15. 

In  the  grid  mode,  the  CTS  system  is  also  coupled  with  the  facility  digital  computer  for 
closed-loop  operation.  However,  in  this  mode,  the  store  model  is  automatically  positioned 
by  the  CTS  at  selected  locations  and  orientations  with  respect  to  the  aircraft  model  (or  in  the 
free  stream)  which  are  pre-programmed  into  the  computer.  Typical  aerodynamic  coefficient 
data  as  obtained  in  Tunnel  4T  are  presented  in  Fig.  16. 

In  many  CTS  grid  mode  applications,  the  store/balance  combination  is  replaced  by  a 
flow  angularity  probe  (Fig.  17).  For  this  case,  local  flow  angles  (Fig.  18)  are  the  data  output 
rather  than  aerodynamic  coefficients.  As  might  be  expected,  most  of  the  experimental  data 
obtained  during  CTS  grid  tests  are  used  in  the  development  of  flow-field  math  models  for 
analytical  trajectory  generation  programs.  However,  at  the  WSRL  (Ref.  11),  experimental 
trajectories  are  generated  using  a  flow  angularity  probe  rather  than  a  store  model /internal 
balance  combination  because  of  the  small  scale  of  aircraft  and  store  models  required  by 
wind  tunnel  size  (typically,  0.02-scale  models  are  used).  The  technique  is  similar  to  other 
experimental  simulations  except  that  each  set  of  store  aerodynamic  coefficient  data  must  be 
calculated  using  flow  angles  measured  at  12  points  along  the  predicted  store  centerline 
location  instead  of  being  directly  measured  with  a  store  model/internal  balance  combination 
placed  in  the  same  position. 
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Advantages  and  disadvantages  of  wind  tunnel  computer  simulations  with  respect  to 
other  trajectory  acquisition  techniques  must  be  considered  in  the  choice  of  techniques.  Since 
this  technique  is  experimental,  many  areas  of  concern  are  hardware  related.  Some  of  the 
advantageous  characteristics  of  the  wind  tunnel  computer  simulation  technique  are 
as  follows: 

1.  Test  results  are  valid  within  the  limits  of  the  simulation.  Herein  lies  probably  the 
single  greatest  advantage  of  this  technique  (with  respect  to  analytical  techniques) 
in  that  measured  store  aerodynamic  data  are  used  to  determine  the  trajectory 
path.  The  only  assumptions  regarding  store  aerodynamic  behavior  in  the  aircraft 
flow  field  are  that  Reynolds  number  and  sting  interference  effects  are  small. 

2.  Test  results  are  available  at  the  completion  of  each  trajectory,  allowing  almost 
real-time  monitoring  of  the  data.  This  markedly  increases  the  flexibility  and 
effectiveness  or  the  test  program  since  subsequent  trajectory  selections  can  be 
based  on  lest  results.  When  unanticipated  trajectory  trends  make  changes  in  test 
program  direction  necessary,  the  required  changes  can  normally  be 
accomplished  quickly  and  easily. 

3.  Several  parameters  can  be  simulated  which  cannot  be  simulated  in  drop  tests. 

These  include  aircraft  acceleration  (pullup/pushover  maneuvers),  aircraft  dive 
angle,  thrust,  and  active  guidance. 

4.  The  number  of  models  required  is  substantially  lower  than  for  drop  tests.  Also, 
design  and  fabrication  of  the  models  is  not  as  difficult  since  only  the  external 
shape  of  the  store  need  be  matched. 

Several  papers  have  been  published  which  discuss  a  number  of  the  problem  areas 
associated  with  wind  tunnel  trajectory  techniques  (e.g.,  Refs.  15  through  17).  Some  of  the 
characteristics  of  this  technique  which  might  be  considered  disadvantageous  are  as  follows: 

1 .  The  use  of  a  wind  tunnel  facility  is  necessary,  and  facility  productivity  must  be 
considered  in  resource  management.  For  example,  approximately  4  to  5 
trajectories  (of  I -sec  duration,  full-scale  time)  can  be  accomplished  per  test  hour 
in  the  AEDC  wind  tunnels.  However,  for  trajectories  of  shorter  duration,  the 
productivity  is  correspondingly  increased. 

2.  Some  store  detail  may  be  lost  at  small  scale.  For  example,  comparison  of 
0.05-and  0.25-scale  models  ol  the  EOGB-II  (Fig.  19)  .indicates  several  items 
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(launch  lugs,  DME  antenna,  umbilical  fitting,  harness  fairing,  side  conduit,  wing 
fences,  and  flap  gaps)  which  are  present  oil  the  larger  model  but  not  on  the 
smaller.  In  some  cases,  lifting  surfaces  (wings,  canards,  strakes)  on  full-scale 
stores  have  airfoil  shapes,  whereas  those  on  small  models  are  normally  flat  plates 
of  the  same  plan  form  with  shaped  leading  and  trailing  edges.  Correct  alignment 
of  lilting  surfaces  is  also  difficult  at  small  scale.  Some  stores  require  external 
shape  modification  in  order  to  accept  a  sting-balance  combination.  This  is  a 
particular  problem  on  stores  with  long,  tapering  afterbodies  such  as  the  MK-84 
LDGP  (Fig.  20). 

Whenever  possible,  it  is  desirable  lo  determine  the  effects  of  reduced  scaling  on 
store  aerodynamics.  When  discrepancies  in  the  free-stream  aerodynamic 
characteristics  of  different  scale  models  are  detected  (for  whatever  reason), 
corrections  can  be  applied  to  the  small-scale  model  data  during  trajectory  testing 
if  so  desired.  For  example,  agreement  among  the  free-stream  data  for  three 
different  scale  slender  bomb  models  (Fig.  21)  is  quite  good  at  Mach  numbers 
greater  than  I,  but  at  subsonic  Mach  numbers,  the  center-of-pressure  location 
on  the  small-scale  models  is  obviously  further  aft.  For  the  subsonic  case,  the 
moment  reference  point  of  the  small-scale  models  could  be  adjusted  to  bring  the 
static  margin  more  in  line  with  the  large-scale  model  data.  However,  in 
implementing  this  adjustment,  it  is  assumed  that  aerodynamic  effects  in  the 
curvilinear  aircraft  flow  held  are  comparable  to  those  in  the  free  stream,  which 
may  or  may  not  be  the  case.  Unpublished  AEDC  data  of  D.  W.  Hill  have  shown 
that  adjustments  to  the  measured  model  data  can  result  in  substantial  changes  in 
trajectory  motion  (Fig.  22). 

3.  For  sting-mounted  models,  the  influence  of  the  sting  on  measured  aerodynamic 
data  is  perhaps  always  in  question,  especially  at  the  large  sting-to-model 
diameter  ratio  (about  0.3  to  0.5)  normally  required  for  structural  integrity  at 
small  scale.  Numerous  reports  have  been  published  documenting  sting  effects  on 
aerodynamic  data  (e.g.,  Refs.  18  through  20).  Cahn  (Ref.  18)  noted  that  the  drag 
characteristics  on  models  with  boattail  afterbodies  were  particularly  susceptible 
to  sting  effects  at  transonic  Mach  numbers  and  that  surface  pressure 
distributions  could  be  altered  as  far  upstream  from  the  base  as  40  percent  of  the 
length  of  the  model.  Changes  in  the  pressure  distribution  on  the  store  afterbody 
would  almost  certainly  affect  the  moment  characteristics  as  well.  The  data  of 
Ref.  18  also  showed  that  sting  interference  on  drag  diminished  rapidly  as  ratio  of 
straight  sting  length  (no  flares  or  bends)  to  base  diameter  increased.  For  a 
straight  sting  length  of  about  six  base  diameters,  interference  effects  were 
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negligible  on  most  boaitail  configurations.  Adherence  to  such  a  criterion  for 
CTS  sting-mounted  models  is  not  always  possible,  it  is  certainly  desirable, 
especially  for  boattailed  stores. 

4.  The  presence  of  the  CTS  mechanism  may  influence  store  aerodynamic  loads. 
Studies  have  been  made  in  Tunnel  4T  to  determine  effects  of  overall  sting  length 
(or  the  proximity  of  the  CTS  head  with  respect  to  the  store  eg)  on  measured 
aerodynamic  data  (Hill  and  Ref.  21).  Typical  data  from  Hill  are  presented  in 
Figs.  23  and  24  which  show  effects  of  CTS  proximity  on  wind  tunnel  Mach 
number  and  measured  aerodynamic  coefficients,  respectively.  On  the  basis  of 
ihese  results,  a  minimum  sting  length  of  24  in.  is  recommended  between  the  CTS 
yaw  center  and  model  center  of  gravity. 

5.  Some  detail  of  the  aircraft  suspension  equipment  may  be  lost  at  small  scale,  such 
as  ventilating  slots  or  sway  braces  on  multiple  racks  or  pylons.  Dix  (Ref.  22)  has 
shown  improved  correlation  of  wind  tunnel  and  flight  test  results  for  MK-83 
carriage  loads  aerodynamics  when  sway  braces  were  present  on  a  triple  ejection 
rack  (TER)  model  (Fig.  25).  Dix  has  also  determined  that  the  presence  of  sway 
braces  could  sometimes  have  a  significant  effect  on  trajectory  motion  (Fig.  26). 
In  view  of  these  results,  it  is  desirable  to  include  as  much  detail  as  possible  on  the 
aircraft  model  suspension  hardware. 

6.  Foreshortened  trajectories  usually  result  for  rapidly  pitching  or  tumbling  stores 
because  of  CTS  hardware  limitations.  The  required  store  position/attitude  in  the 
aircraft  flow  field  may  cause  the  store  model  support  sting  or  CTS  to  contact  the 
aircraft  or  its  support  sting.  Also,  to  set  large  trajectory  attitudes  may  require  a 
CTS  position  outside  of  the  physical  travel  limits  even  though  the  store  may  still 
be  in  close  proximity  to  the  aircraft. 

7.  Store  rolling  capability  is  sacrificed  when  many  offset  and  bent  sting-support 
combinations  are  used  in  order  to  provide  adequate  clearance  between  the  CTS- 
and  aircraft-related  hardware.  How  the  loss  of  this  degree  of  freedom  is  dealt 
with  in  the  Tunnel  4T  and  VKF  trajectory  programs  is  contained  in  Volumes  II 
and  Ill  of  this  series,  respectively. 

8.  The  spacing  between  the  store  model  and  suspension  equipment  at  trajectory 
initiation  can  affect  trajectory  results.  For  a  vibration-free  environment,  a  gap 
between  the  model  surface  and  suspension  equipment  of  nominally  0.070  in., 
model  scale  (approximately  1.4  in.,  full  scale)  exists  when  the  store  model  is 
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located  at  the  carriage  position.  This  roughly  approximates  full-scale  values. 
However,  because  of  vibration  of  the  store  or  aircraft  models,  it  is  sometimes 
necessary  to  initiate  trajectories  at  a  greater  distance  from  the  carriage  location 
to  prevent  premature  collisions.  Since  aerodynamic  gradients  are  usually  large 
near  carriage,  nonrepresentative  trajectories  could  result.  It  is  generally 
recommended  that  the  gap  between  store  and  aircraft  suspension  equipment  be 
no  greater  than  0.150  in.,  model  scale  (approximately  3  in..  Full  scale)  at 
trajectory  initiation. 

9.  Correcting  the  model  attitude  for  angle  of  attack  and  sideslip  angle  induced  by 
the  motion  of  the  full-scale  store  may  produce  changes  in  the  interference 
aerodynamics.  Because  of  the  aerodynamic  angles  induced  by  lateral  and 
vertical  velocities  of  the  store  eg,  the  model  cannot  be  positioned  at  the  correct 
attitude  with  respect  to  both  the  relative  wind  vector  and  the  aircraft.  For  the 
Tunnel  4T  trajectory  program,  the  angular  position  with  respect  to  the  relative 
wind  vector  was  deemed  more  important,  and  the  store  model  is  positioned 
accordingly.  For  the  VKF  trajectory  program,  the  model-induced  angle 
corrections  are  recommended,  but  can  be  deleted  as  a  program  option. 
Fortunately,  store  velocities  in  the  aircraft  near  flow  field  are  usually  not  large; 
thus  the  small  corrections  for  induced  angles  may  result  in  only  minimal  errors 
in  interference  aerodynamics.  At  greater  distances  from  carriage,  store 
velocities  (hence,  induced  angles)  can  be  large,  but  substantial  changes  in  store 
attitude  when  it  is  located  in  the  aircraft  far  flow  field  may  produce  only 
negligible  errors  in  the  interference  aerodynamics.  The  induced  angle 
corrections  for  the  Tunnel  4T  and  VKF  trajectory  programs  are  contained  in 
Volumes  II  and  III  of  this  series,  respectively. 

10.  The  aerodynamic  damping  characteristics  of  most  stores  are  not  well  defined 
even  in  the  free  stream  and  may  well  be  sensitive  to  position  and  attitude  in  the 
highly  nonuniform  flow  field  surrounding  a  typical  aircraft.  Because  of  the  lack 
of  knowledge  of  damping  values,  aerodynamic  damping  derivative  coefficients 
are  input  as  constant  values  in  most  trajectory  programs.  However,  this 
disadvantage  is  not  only  related  to  experimental  trajectory  programs  but  is  a 
factor  in  computer  simulations  in  general.  Marshall  and  Summers  (Ref.  23) 
determined  that  a  ±  50-percent  variation  in  constant  values  of  damping 
produced  little  change  in  trajectory  motion  for  low  amplitude  angular 
excursions.  However,  Hill  has  shown  that  variations  in  aerodynamic  damping 
can  produce  large  differences  in  trajectory  results  for  high  amplitude  angular 
motion  (Fig.  27). 
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3.4  ANALYTICAL  APPLICATIONS 

As  previously  mentioned,  the  basic  difference  between  analytical  and  experimental 
trajectory  programs  is  in  the  method  of  definition  of  the  aerodynamic  coefficient  data.  For 
the  analytical  technique,  aerodynamic  characteristics  of  the  released  store  are  extracted  from 
a  flow-field  math  model.  The  basic  options  available  in  the  design  of  such  a  math  model  are 
described  in  Fig.  1.  In  general,  free-stream  aerodynamics  are  combined  with  interference 
aerodynamics  induced  by  the  aircraft  flow  field  to  define  the  store  total  aerodynamic 
characteristics,  which  are  then  input  into  trajectory  calculations.  Either  theoretical 
predictions  or  an  empirical  data  base  may  be  used  to  define  free-stream  and  interference 
aerodynamics.  However,  to  make  accurate  predictions,  the  flow-field  model  needs  to 
consider  the  aerodynamic  shape  of  the  aircraft  and  store,  aircraft  angle  of  attack  and 
sideslip  angle,  aircraft  external  store  Ipading,  free-stream  Mach  number,  and  the  position 
and  attitude  of  the  store  relative  to  the  aircraft. 

The  formulation  of  a  flow-field  math  model  from  an  empirical  data  base  is  usually 
accomplished  along  the  following  guidelines.  The  free-stream  aerodynamic  characteristics 
are  obtained  from  test  data  for  the  largest  scale  model  available,  since  these  data  are  more 
likely  to  be  representative  of  full-scale  store  aerodynamics.  Interference  aerodynamics  may 
be  determined  by  the  interference  coefficient  method  or  the  flow  angularity  method. 

In  the  interference  coefficient  method,  a  matrix  of  aerodynamic  data  is  measured  on  a 
small-scale  model  at  selected  positions  and  attitudes  in  the  aircraft  flow  field  which  are 
expected  to  encompass  the  subsequent  trajectory  path.  Then  free-stream  aerodynamics  (for 
the  same  store  model  at  the  same  attitude)  are  subtracted  from  the  flow-field  aerodynamics 
to  define  the  interference  coefficients.  In  the  process  it  is  hoped  that  any  effects  of  loss  of 
detail  at  small  scale  pn  store  aerodynamics  will  be  minimized. 

In  the  flow  angularity  method,  flow  angles  are  obtained  with  a  flow  angle  probe  at 
selected  positions  in  the  aircraft  flow  field.  Using  the  flow  angles  relative  to  the  store  body 
and  lifting  surfaces,  one  can  calculate  interference  aerodynamic  coefficients. 

Once  the  interference  coefficients  have  been  defined,  total  store  aerodynamics  used  in 
the  trajectory  simulation  are  determined  at  each  point  along  the  trajectory  path  by  adding 
the  free-stream  (function  of  attitude  with  respect  to  the  relative  wind)  and  interference 
(function  of  position  and  attitude  with  respect  to  the  aircraft)  contributions.  (Although  this 
description  is  intended  primarily  for  empirical  data,  a  math  model  using  theoretical 
estimates  for  the  aerodynamic  coefficients  or  aerodynamic  flow  angles  could  be  assembled 
in  the  same  manner.) 
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A  number  of  analytical  trajectory  programs  have  been  developed  for  separation  analysis. 
Nielsen  Engineering  and  Research,  Inc.  (Refs.  24  through  26)  used  analytic  computations 
such  as  slender  body  theory,  the  vortex-lattice  method,  and  three-dimensional  source 
distributions  in  the  formulation  of  their  flow-field  model.  Rockwell  International  (Ref.  27) 
and  the  Air  Force  Armament  Laboratory  (ADTC/DLJC)  (Ref.  28)  have  developed  flow- 
field  models  based  on  the  flow  angularity  method.  The  ADTC/DLJC  (Ref.  29)  and  Sandia 
Laboratories  (Ref.  30)  have  developed  flow-field  math  models  based  on  the  interference 
coefficient  method.  At  the  AEDC,  most  analytical  trajectory  programs  use  the  interference 
coefficient  method.  Many  other  facilities  have  working  analytical  trajectory  programs;  those 
mentioned  were  selected  simply  to  illustrate  the  diversity  in  the  formulation  of  flow-field 
math  models. 

Advantages  and  disadvantages  of  the  analytical  computer  simulations  with  respect  to 
other  trajectory  acquisition  techniques  must  be  considered  in  their  development  and  use. 
Some  of  the  characteristics  of  the  analytical  computer  simulation  technique  are  as  follows: 

1.  When  compared  to  experimental  techniques,  the  cost  per  trajectory  for 
analytical  simulations  can  be  significantly  lower.  Obviously,  considerable 
engineering  and  programming  expense  is  involved  in  the  development  of  an 
analytical  program;  however,  once  a  program  is  developed,  large  numbers  of 
trajectories  can  be  calculated  at  a  relatively  small  cost  per  trajectory. 

2.  Analytical  simulations  are  quite  flexible.  Within  the  limitations  of  the  flow-field 
model,  the  influence  of  a  number  of  parameters  may  be  investigated.  In  fact, 
parametric  evaluations  are  ideally  suited  to  this  technique. 

3.  As  for  the  wind  tunnel  computer  technique,  several  parameters  such  as  aircraft 
acceleration,  aircraft  dive  angle,  thrust,  and  active  guidance  can  be  simulated 
which  cannot  be  simulated  in  drop  tests. 

4.  Trajectory  results  are  immediately  available.  No  delays  are  inherent  in  this 
technique  which  might  affect  design  decisions,  flight  clearance  status,  or 
subsequent  analysis  effort. 

5.  Analytical  simulations  can  be  accomplished  which  require  no  empirical  data. 
However,  as  evidenced  by  the  Nielsen  program  (Refs.  24  through  26), 
development  of  such  simulations  may  require  years  of  work.  A  complete 
description  of  nonuniform  flow  fields  around  complex  configurations  is  not 
amenable  to  theoretical  techniques  at  this  time. 
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6.  The  flow  angularity  method  is  perhaps  slightly  more  flexible  than  the 
interference  coefficient  method.  Once  the  aerodynamic  flow  angles  in  an  aircraft 
flow  field  are  defined,  trajectories  for  different  shaped  stores  can  be  determined 
with  minor  adjustments  to  the  flow-field  calculations,  and  no  additional 
empirical  data  are  required  (assuming  the  aircraft  loading  configuration  is  not 
altered).  However,  of  all  of  the  analytical  flow-field  models,  those  using  the 
incremental  coefficient  method  are  probably  the  easiest  to  generate. 

7.  The  adaptation  of  this  technique  for  the  extension  of  experimentally  obtained 
trajectories  which  were  prematurely  terminated  because  of  CTS  hardware 
limitations  is  very  feasible  (Fig.  28).  Development  of  this  capability  at  the  AEDC 
is  currently  in  progress,  and  documentation  should  be  completed  in  the  near 
future.  In  many  cases,  extensions  of  this  nature  could  provide  the,  basis  for 
additional  testing  or  for  release  certification. 

8.  It  is  extremely  difficult  to  develop  a  comprehensive  flow-field  model  which 
accurately  accounts  for  all  variables.  For  -this  reason,  the  reliability  of 
trajectories  generated  by  analytical  programs  is  not  yet  on  a  par  with  wind  tunnel 
techniques  and  may  be  suspect  unless  demonstrated  otherwise.  It  should  be 
noted  that  good  results  have  been  obtained  with  this  technique  in  many  cases  for 
limited  flow-field  applications  (for  example,  see  Figs.  22  and  27). 

4.0  CONCLUDING  REMARKS 

As  the  sophistication  of  modern  weaponry  is  increased  and  release  certification  is 
required  at  more  demanding  flight  conditions,  the  complexity  of  store  separation  analysis 
programs  must  necessarily  increase.  Thus,  more  thorough  planning  will  be  required  to 
effectively  util^e  the  available  prediction  techniques  in  the  attainment  of  overall  program 
objectives.  In  the  formulation  of  a  test  program,  available  resources  (manpower  for  design, 
model  fabrication,  programming,  and  analysis;  access  to  test  and  computer  facilities; 
dollars;  etc.)  must  be  balanced  against  the  advantages  and  disadvantages  of  each  trajectory 
technique  to  arrive  at  the  combination  that  will  accomplish  the  objectives  of  the  test 
program  in  the  most  timely  and  efficient  manner.  Obviously,  this  combination  is  almost 
certain  to  be  unique  for  each  specific  application.  In  the  discussion  of  the  trajectory 
techniques  in  Sections  2  and  3,  an  attempt  was  made  to  present  the  benefits  and  limitations 
of  each  in  such  a  way  as  to  demonstrate  their  truly  complementary  nature,  and,  in  so  doing, 
to  provide  information  which  may  be  used  for  guidance  in  the  formulation  of  future  store 
separation  lest  and  evaluation  programs.  Specifically,  use  of  the  analytical  simulation 
technique  is  encouraged  whenever  possible  to  aid  in  structuring  efficient  and  cost  effective 
wind  tunnel  test  programs  so  that  maximum  benefit  may  be  realized  from  the  empirical  data. 
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Figure  1.  Diagram  of  trajectory  generation  options. 
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Figure  2.  Three-DOF  coordinate  notation. 
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To  augment  the  ejector  force  for  insufficient  gravity  up  to 
the  end  of  stroke : 
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Figure  3.  Summary  of  dynamic  similarity  principles. 
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a.  Incremental  flow  angle 
Figure  5.  Aircraft  flow  field  for  stable  store. 
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b.  Incremental  pitching-moment  coefficient 
Figure  5.  Concluded. 
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a.  Transonic  aerodynamics 
Figure  6.  Froude  scaling  with  stable  store. 
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a.  AFe!  =  0 

Figure  8.  Light  Mach  scaling  with  stable  store. 


-  I  0 


b.  AF^  =  m'  (Ag') 
Figure  8.  Continued. 
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c.  AFe!  =  3m'  (Ag') 
Figure  8.  Concluded. 
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YES 


Figure  10.  Simplified  block  diagram  of  a  trajectory 
generation  computer  program. 
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Figure  13.  Isometric  drawing  of  a  typical  store  separation  installation  in 

Tunnel  4T  and  a  block  diagram  of  the  computer  control  loop. 


AEDC-TR-79-1 


Figure  14.  Isometric  drawing  of  a  typical  store  separation 
installation  in  Tunnels  B  and  C. 
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at, deg 


a,  deg 


a,  Carriage  position  characteristics 
Figure  16.  Typical  aerodynamic  coefficient  data  obtained 
on  a  slender-finned  bomb  in  the  flow  field  of 
a  supersonic  fighter  aircraft. 
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Figure  18.  Flow-field  measurements  beneath  a  supersonic 
fighter  aircraft. 
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Figure  21.  Variation  of  the  zero-lift  longitudinal  stability 
derivatives  and  axial-force  coefficients  of  a 
slender-finned  bomb  with  Mach  number. 
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Figure  22.  Effects  of  coefficient  corrections  on  a  typical 

trajectory  of  a  slender-finned  bomb. 
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WITHOUT  SWAY  BRACES 
WITH  SWAY  BRACES 


Ref.  22 


Figure  26.  Effect  of  sway  braces  on  separation  trajectories 

of  the  MK-83  low-drag  bomb  from  the  F-4C  Aircraft, 
left  inboard  pylon,  TER  station  1,  level  flight, 

1,200-lb  ejector  force. 
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The  angular  momentum  of  a  rotating  rigid  body  about  a  point  other  than  the  mass  center 
may  be  determined  as  follows1: 


Oe(XeYeZe)  is  a  space-fixed  axis  system. 


Xg  is  parallel  to  Xg,  Yg  is  parallel  to  Yg,  and  Zg  is  parallel  to  Zg  at  all  times. 


is  the  velocity  of  point  0  (with  respect  to  Og). 


-> 

6J 


is  the  angular  velocity  about  point  0. 


Vi  =  Vo  +  ^X  ri 


+  r 


The  angular  momentum  of  the  body  about  point  O  is: 

H  =  •  x  m iV- 

O  111 


’Etkin,  Bernard.  Dynamics  of  Flight:  Stability  and  Control.  John  Wiley  and  Sons,  Inc.,  New 
York,  1959,  pp.  94-99,  442-443. 
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Substituting  for  r;  and  v;,  and  expanding  (noting  that  £rmj  =  0)  gives: 

mr  xv  +  mr  x  (&>  x  r  )  +  2,  r  x  m-(oj  x  r) 

O  O  o  U  1 

However,  the  last  term  in  the  above  relation  is  defined  as  the  angular  momentum  about  the 
mass  center  (H),  such  that 

Hq  =  mrQ  x  vQ  +  mr0  x  (co  x  rQ)  +  H 
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APPENDIX  B 

GENERALIZED  MOMENT  EQUATION  AND  ANGULAR  ACCELERATIONS 


For  a  nonrotating  XpYuZuaxis  system  (as  shown  in  Appendix  A),  the  moment  equation 
is: 


£  M 

O 


£(JU 

dt 


However,  by  letting  the  XbYbZb  axis  system  rotate  (at  angular  velocity  0),  one  obtains 
certain  advantages,  and  the  moment  equation  becomes: 


->  8  (H  ) 

^0~^xHo 


Substituting  for  H0  (from  Appendix  A)  and  expanding  gives 

S  (H)  S 


1  M°  =  TT  +  sT  {m"r°  x  ''°)  +  sT  C"1"0  x  x  '->] 

+  n  x  H  +  n  X  (mro  X  vQ)  +  Ox  (mr0  x  (<u  x  7Q)  J 


For  most  store  separation  applications,  the  following  criteria  are  valid: 
1.  The  angular  momentum  about  the  mass  center  is  H  =  Iu> 


[xx 

_IXY 

_IXZ 

where 

I  = 

-IXY 

IYY 

_IYZ 

_IXZ 

-IYZ 

lZZ 

2.  Fixing  the  rotating  axis  system  with  respect  to  the  body  gives 


0  =  <y  =  ip  +  j  q  +  kr 

MO)  S(<7)  t.  1.  7. 

-  =  -  =  ip  +  iq  +  kr 

St  St 

i  =  o 
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3.  Aligning  the  XB-ZB  plane  of  the  axis  system  in  the  body  plane  of  symmetry  gives 

*XY  =  *YZ  =  0 

4.  The  distance  from  the  rotation  center  to  the  mass  center  (?0)  is  nonzero  only  during 
staged  release,  but  in  any  case,  r0  is  not  time  dependent. 

r  =  T  X  +  T  Y  +  k  Z 

O  O  J  o  o 


S(r  ) 


-=  0 


5.  When  r0  st  0,  the  velocity  and  acceleration  of  the  rotation  center  are 

'o,Z 


vo  =  1  vo,X+  i  vo, Y  +  k 


5(v  ) 


Cs  —  1  3.  y  1  3  _  y  t  It  fct  ry 

Ot  °,  A  Jo,i  o,  Z, 

6.  If  rt,  sfc  0,  the  inertia  about  the  rotation  center  may  be  defined  as 

^  +  Zjj) 


+  t  ar 


*xx  =  *xx 


o  o 


*XY  “  *XY  +  mXoYo 


‘XZ 


IY  7  +  mX  Z 

A  Z  0  0 


M- 


11 

+  m(X2  + 

YZ  =  *YZ 

+  mY  Z 

o  o 

zz  =  *zz 

+  m(X2  + 

Vo  ) 

rxx 

_IXY 

-IXZ 

~Txy 

iyy 

_IYZ 

_IXZ 

-IYZ 

*zz 

70 


AEDC-TR-79-1 


7.  The  scalar  components  of  M,„  H,  and  5(H)/<5t  are  defined  as 

-*  -> 

M  =  i  M  Y  +  j  M  v  +  kM  7 

O  0,A  J  0,1  o,Z 


Hx 

p 

hy 

-[T] 

q 

Hz 

r 

«X 

P 

hy 

=  [T] 

q 

Hz 

r 

By  applying  the  preceding  criteria,  one  can  determine  the  scalar  components  of  the  moment 
equation  to  be: 

2Mo,X  =  HX  +  qHz  -  rHy  +  mXo(qvo  Y  +  rvo>z) 

+  “Vvz-^o.x)  -  mZoK,Y  +  rvo,X^ 

SMq  Y  =  Hy  +  rHx-pHz -mXo(ao  Z  +  pvo  Y) 

+  mYo(Pvo,X  +  rvo,Z^  +  mZo(ao,X  -  rv0> 

YMo,Z  =  HZ  +  PHY~qHX  +  mXo(ao,Y  -  Pvo,Z) 

“mYo(ao,X  +  V'o.Z*  +  mZo(Pvo,X  +  qvo,Y> 

The  moments  about  the  Xu,  Yu,  and  Zb  axes  at  the  rotation  center  are  (see  Fig.  B-l): 

SMo,X  =  Mx  +  FzY0-FyZo  +  R£ 

2MoY  =  My  +  FxZo-FzXo  +  Rm 

2Moz  =  Mz  +  FyXo-FxYo  +  Rn 
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Substituting,  then: 


11  x  =  Mx  +  fzYo~fyZo  +  R£~qHz  +  rHY 

-mXo(qvb)Y  +  rv0)Z)-niYoK,Z -qvo,x)  +  mZo(ao,Y  +  rvo,x) 

HY  =  My  +  FxZo-FZXo  +  Rm-  rHX  +  PHZ 

+  mXo(ao,Z  +  Pvo,Y)-mYo(Pvo,X  +  rvo,Z}  -mZo(ao,X-rvo,Y) 

Hz  =  Mz  +  FYXo-FXYo  +  Rn~PHY  +  ?HX 

-mXo(ao,Y  ~Pvo,z)  +  mYo(ao,X  +  qvo,z)  -  mZo(PVo,X  +  cIvo)Y) 

The  angular  accelerations  are  then: 


p 

Rx 

r— ■ i  —  i 

q 

=  H 

Hy 

r 

where  [I]-1  is  the  inverse  of  [!]. 
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Figure  B-1.  Positive  directions  of  full-scale  forces  and  moments. 
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APPENDIX  l 

GENERALIZED  FORCE  EQUATION  AND  LINEAR  ACCELERATIONS 


The  generalized  force  equation  is: 


0^(X[;Y l;Zh)  is  fixed  frame  of  reference. 

0(XiiY|jZH)  is  second  frame  of  reference  moving  with  a  linear  velocity  (vG)  at  the  origin  and 
at  angular  velocity  (u>)  with  respect  to  the  fixed  frame. 

The  velocity  at  the  mass  center  (v)  may  be  determined  as: 

-»  d  (r  ''  ) 

v  =  - ; - 

d.t 

d  (r  " )  d  (r  ) 

o  o 

—  — — —  -f-  “ 

dt  dt 

->  3  /->  >.  -> 

_  V°  +  8t  r°  +  “  x  r° 
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By  defining: 


v  =  i  u  +  j  v  +  kw 


%  =  ivo,X  +  ivo,Y  +  kvo,Z 


r„  =  ?X  +  JY  +  kZ 

o  o  J  o  o 


(o  =  ip  +  jq  +  kr 


The  scalar  components  of  velocity  are: 


u  =  v  x  +  qZo-rYc 


v  =  vo,  Y  +  rXo  ~  P^o 


w  =  yoZ  +  pYo-qX0 

The  acceleration  of  the  mass  center  may  be  determined  to  be: 


+  a>  x  v 


Substituting  for  v  and  expanding  gives: 


fi(v  )  S2(v  )  8(2)  x^  ) 

n  u  u 


+  0)  x  vo  +  0)  X  — - —  +  CO  x  (O)  x  rQ) 
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Then  by  definition  one  arrives  at: 

?  ■?  p 

3=1  3-j^  +  J3y  + 

S(vQ)  ->  _»  -> 

-jr-  =  iao,X  +  jao,  Y  +  kao,  z 

S2(r) 

- 2-  =  0 

St2 

5(<m)  ■?.  “?•  j*- 

— —  =  ip  +  jq  +  kr 
St 

The  scalar  components  of  acceleration  are  then: 

aX  =  ao,X  +  1vo,Z  _  rvo, Y 

-XQ(q2  +  r2)  +  Y0(pq-r)  +  Z0(q  +  pr) 

aY  =  ao,Y  +  rvo,X  “  Pvo, Z 

+  Xo(pq  +  r)  -  Yo(p2  +  r2)  +  zQ(qr  —  p) 

aZ  =  ao,Z  +  Pvo,Y  -  <Ivo,X 

+  XQ(pr  -  q)  +  YQ(qr  +  p)  -  ZQ(p2  +  q2) 
-> 

The  scalar  components  of  IT  are  (see  Fig.  B-l): 

IF  =  i(Fx  +  Rx)  +  j  (F  y  +  Ry)  +  k(Fz  +  Rz) 

Therefore  the  scalar  force  equations  are: 

Fx  +  Ry  =  max 

Fy  +  By  —  nt a y 

FZ  +  P'Z  =  maz 
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u  =  ao  x,  V  =  ao  Y,  vV  =  ao  Z 

The  linear  accelerations  of  the  mass  center  arc  then: 

u  =  Fx/m  +  rv  —  qw 
'm  =  F  y/m  +  pw  -  ru 
w  =  F^/m  +  qu  -  pv 

For  a  staged  release  (r0  =£  0),  the  terms  u,  v,  w  do  not  appear  in  the  acceleration  equations. 
However,  by  substituting  one  obtains: 


a^  =  u  +  qw 


ay  =  v  +  ru  —  pw 

a  2  =  w  +  pv  —  qu 

vo,X=  u  -qZo  +  rYo 
vo,Y  =  v  —  rXQ  +  pZo 

vo,z  =  w-pyo  +  qxo 

Then,  the  linear  accelerations  are: 

u  =  a  Y  +  qZ  —  rY 
o,A  *  o  o 

v  =  a  v  +  rX  —  pZ 
o,  i  o  r  o 

W  =  ao,Z  +  PYo-^Xo 

Note  that  for  a  staged  release,  the  linear  accelerations  of  the  mass  center  are  influenced 
by  the  applied  forces  only  through  coupling  with  the  angular  accelerations. 
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APPENDIX  1) 

FULL-SCALE  FORCE  AND  MOMENT  EQUATIONS 

For  a  free-falling  store,  components  of  the  full-scale  forces  and  moments  acting  through 
or  about  the  store  eg  arc  described  as  follows: 

Forces 

FX  =  WX-(l~A(CA,t  +  CA,o)  +  FT,X  +  FE,X 

FY  =  ^  Y  +  1x>A^Y  +  CY,o)  +  FT,Y  +  FE,Y 

FZ  =  WZ_<?°oA(CN  +  CN,0)  +  FT,Z  +  FE,Z 

Moments 

MX  =  ^ooA^3  [C£  +  C£,0  +  C(!p(^pV2Ur)]  +  mt,x  +  me,x 

MY  =  <UAei[Cm  +  Cm,o  +  Cm(i(q£l/2UR)]  +  MT,Y  +  ME,Y 

MZ  =  ^[c.  +  Cn  o  +  (^r£2/2UR^J  +  Mt  z  +  Me  z 

Additional  terms  required  to  define  the  full-scale  moments  and  forces  which  occur  on  staged 
release  are  given  in  Appendixes  B  and  C,  respectively. 
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A 

■» 

a 


a0,x,ao,Y 

a0,z 

ax.^Yi^z 


aZ,p 

Ca,i.Cn,Cy 


Ca,oi  Cn  ,0* 
Cy,o 

C|\Cm,Cn 


Cn  rv 


Ci’p.Cjnq.Cn,. 


ACn 


Cmc 

Cn,. 


NOMENCLATURE 

Store  full-scale  reference  area,  ft2 

Total  acceleration  of  the  store  mass  center,  ft/sec2 

Accelerations  of  the  store  rotation  center  (if  not  coincident  with  the 
mass  center)  in  the  XB,  YB,  and  ZB  directions,  respectively,  ft/sec2 

Accelerations  of  the  store  mass  center  in  the  XB,  YB,  and  ZB 
directions,  respectively,  ft/sec2 

Acceleration  of  the  aircraft  in  the  Zf  direction,  ft/sec2  (see  Fig.  12) 

Store  measured  (or  predicted)  axial-force,  normal-force,  and  side- 
force  coefficients,  positive  in  the  negative  XB,  negative  ZB,  and 
positive  YB  directions,  respectively 

External  input  axial-force,  normal-force,  and  side-force  coefficients, 
respectively 

Store  measured  (or  predicted)  rolling-moment,  pitching-moment, 
and  yawing-moment  coefficients,  respectively;  the  positive  vectors 
are  coincident  with  the  positive  XB,  YB,  and  ZB  directions 

External  input  rolling-moment,  pitching-moment,  and  yawing- 
moment  coefficients,  respectively 

Store  roll-damping,  pitch-damping,  and  yaw-damping  derivatives, 
respectively,  per  radian 

Incremental  value  of  pitching-moment  coefficient  used  in  the  Ref.  8 
flow-field  math  model 

Pitching-moment  coefficient  derivative,  per  radian  or  per  degree 
Normal-force  coefficient  derivative,  per  radian  or  per  degree 
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dCm/ <JCn 


Static  margin,  reference  lengths  forward  of  the  moment  reference 
center 


F 

Fei 

AFe, 


Fe,x.Fe,y» 

Fe,z 


Total  force  acting  on  the  store,  lb 
Total  ejector  force  acting  on  the  store,  lb 

Model  ejector  force  augmentation  to  compensate  for  gravity 
deficiency,  lb 

Components  of  the  ejector  force  acting  on  the  store  in  the  positive 
XB,  Yb,  and  Zb  directions,  respectively,  lb 


Ft,x.Ft,y» 


T,Z 


Fx,Fy,Fz 


Components  of  the  thrust  force  acting  on  the  store  in  the  positive  XB, 
Yb,  and  Zb  directions,  respectively,  lb 

Components  of  total  force  acting  on  a  free-falling  store  in  the 
positive  Xb,  Yb,  and  Zb  directions,  respectively,  lb 


Ag' 

H 


Acceleration  of  gravity,  ft/sec2 

Model  gravitational  deficiency,  ft/sec2 

Angular  momentum  of  a  body  rotating  about  its  mass  center, 
ft-lb-sec 


H0 


Angular  momentum  of  a  body  rotating  about  a  point  other  than  its 
mass  center,  ft-lb-sec 


Hx,Hy,Hz 


Scalar  components  of  the  angular  momentum  (H)  about  the  XB, 
Yb,  and  ZB  axes,  respectively,  ft-lb-sec 


h 

I 


Simulated  pressure  altitude,  ft 

Notation  for  the  inertia  matrix  of  a  store  which  is  calculated  about 
the  mass  center  (see  Appendix  B) 


Notation  for  the  inertia  matrix  of  a  store  which  is  calculated  about  a 
point  (rotation  center)  other  than  the  mass  center  (see  Appendix  B) 
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IxxJyyJzz 


IxxJyyJzz 


Ixy.IxzJyz 


IxyJxzJyz 


*  t 

i.J 


Full-scale  moments  of  inertia  about  the  store  Xb,  Yb,  and  Zb  axes, 
respectively,  slug-ft2,  and  referenced  to  the  store  mass  center 

Full-scale  moments  of  inertia  about  the  store  Xb,  Yb,  and  Zb 
directions,  respectively,  slug-ft2,  but  referenced  to  a  point  (rotation 
center)  other  than  the  mass  center 

Full-scale  products  of  inertia  in  the  store  XB-YB,  Xb-Zb,  and  YB-ZB 
planes,  respectively,  slug-ft2,  and  referenced  to  the  store  mass  center 

Full-scale  products  of  inertia  in  the  store  Xb-Yb, Xb-Zb,  and  YB- 
ZB  planes,  respectively,  slug-ft2,  but  referenced  to  a  point 
(rotation  center)  other  than  the  mass  center 

Unit  vectors  along  the  store  XB,  Yb,  and  Zb  directions,  respectively 


?p  Probe  extension  length,  in.  (see  Fig.  23) 

fi,  t2,  Store  full-scale  reference  dimensions  for  pitching-moment,  yawing- 

moment,  and  rolling-moment  coefficients,  respectively,  ft 

Me  x,Me(y,  Components  of  ejector  moment  about  the  Xb,  Yb,  and  Zb  axes, 

MEjz  respectively,  ft-lb 

-> 

M0  Total  moment  acting  on  the  store,  ft-lb 


M0ix,M0jY,  Scalar  components  of  the  total  moment  acting  on  the  store  about  the 

M0>z  XB,  Yb,  and  ZB  axes,  respectively,  ft-lb 


Mt,x>Mt,y>  Components  of  the  thrust  moment  acting  on  the  store  about  the  Xb, 

Mt>z  Yb,  and  Zb  axes,  respectively,  ft-lb 

Mx,My,Mz  Components  of  the  total  moment  acting  on  a  free-falling  body  about 

the  XB,  Yb,  and  ZB  axes,  respectively,  ft-lb 

M*,  Free-stream  Mach  number 

AM  a,  Free-stream  Mach  number  increment 
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m 

mi 


Nz 


P.q.r 


qP 


Qo= 

Rf>Rm>Rn 


RP 


Rx,  Ry>  Rz 


-> 

r 


-► 

r 


n 


-» 


Tt 


T 


oo 


Store  mass,  slugs 

Store  elemental  mass,  slugs 

Aircraft  “g”  loading  factor 

Store  angular  velocity  about  the  XB,YB,  and  Zb  axes,  respectively;  the 
positive  vectors  are  coincident  with  the  XB,  YB,  and  ZB  axes,  rad/sec 

Pitch  rate  of  the  aircraft  during  a  pullup/pushover  maneuver, 
rad/sec  (see  Fig.  12) 

Free-stream  dynamic  pressure,  psf 

Full-scale  body-axis  restraining  moments  about  the  pivot  (rotation 
center);  the  positive  vectors  are  coincident  with  the  positive  XB,  YB, 
and  ZB  directions,  respectively,  ft-lb 

Effective  rotation  arm  of  the  aircraft  during  a  pullup/pushover 
maneuver,  ft  (see  Fig.  12) 

Full-scale  body-axis  pivot  (rotation  center)  restraining  forces, 
positive  in  the  positive  XB,  YB,  and  ZB,  directions,  respectively,  lb 

Distance  from  the  mass  center  to  an  elemental  mass,  ft 

Distance  from  the  fixed  frame  of  reference  to  the  mass  center,  ft 

Distance  from  the  rotation  center  to  an  elemental  mass,  ft 

Distance  from  the  rotation  center  to  the  mass  center,  ft 

Distance  from  the  fixed  frame  of  reference  to  the  rotation  center,  ft 

Free-stream  total  temperature,  °R 

Free-stream  static  temperature,  °R 
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Trajectory  time  from  the  instant  of  store  release  from  the 
aircraft,  sec 


Time  of  action  of  the  model  ejector,  sec 


UA 

Ur 


Velocity  of  the  aircraft  at  the  simulated  altitude,  ft/sec 

Total  velocity  of  the  full-scale  store  with  respect  to  a  space-fixed 
point,  ft/sec 


u,v,w 


Vn 


Velocities  of  the  full-scale  store  relative  to  the  origin  of  the  flight-axis 
system  in  the  positive  XB,  YB,  and  ZB  directions,  respectively,  ft/sec 

Free-stream  velocity,  ft/sec 


Total  velocity  of  the  store  mass  center,  ft/sec 

Velocity  of  the  elemental  mass  with  respect  to  a  space-fixed  point, 
ft/sec 

Velocity  of  the  rotation  center  with  respect  to  a  spaced-fixed  point, 
ft/sec 


V0,X.V0,Y, 

v0,z 


Scalar  components  of  rotation  center  velocity  (v0),  positive  in  the 
positive  XB,  YB,  and  ZB  directions,  respectively,  ft/sec 


Wt 

WX,WY,WZ 


Full-scale  store  weight,  lb 

Components  of  the  full-scale  store  weight  acting  in  the  positive  XB, 
Yb,  and  ZB  directions,  respectively,  lb  (including  effects  of  simulated 
dive  or  bank  angles) 


X,Y,Z 


Separation  distance  of  the  store  eg  from  the  flight-axis  system  origin 
in  the  positive  XF,  YF,  and  ZF  directions,  respectively,  ft,  full  scale 


X, 


eg 


Axial  distance  from  the  store  nose  to  the  eg  location,  ft,  full  scale 


X 


El 


Distance  from  the  mass  center  to  the  ejector  line  of  action  in  the  XB 
direction,  ft,  positive  forward 
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X[,Yi,Z[  Separation  of  the  store  eg  from  the  inertial-axis  system  origin  in  the 

positive  Xj,  Y1;  and  Z]  directions,  respectively,  ft,  full  scale 

X0,Y0,Z0  Distances  from  the  pivot  (rotation  center)  to  the  store  center  of 

gravity  along  the  Xg,  Yb,  and  Zb  axes,  respectively,  positive  in  the 
positive  XB,  YB,  and  ZB  directions,  ft,  full  scale 

Xp,YP,Zp  Separation  distance  of  the  store  eg  from  the  pylon-axis  system  origin 

in  the  positive  Xp,  Yp,  and  Zp  directions,  respectively,  ft,  full  scale 

AXm,cg  Axial  distance  from  the  store  eg  to  the  pitching-moment  reference 

center,  positive  in  the  positive  XB  direction,  ft,  full  scale 

AZ'  Calculated  adjustment  to  the  observed  store-model  vertical 

displacement  to  correct  for  gravity  deficiency  and  ejector-force 
augmentation  effects,  ft 

a  Aircraft  model  angle  of  attack  relative  to  the  free-stream  velocity 

vector,  deg 

as  Store  model  angle  of  attack,  deg 

axY  Sidewash  angle  with  respect  to  the  store  longitudinal  axis,  positive 

inboard  as  seen  by  the  pilot  for  aircraft  left-wing  data,  deg 

«xz  Upwash  angle  with  respect  to  the  store  longitudinal  axis,  positive  up 

as  seen  by  the  pilot,  deg 

Aas  Incremental  angle  of  attack  on  the  store  used  in  the  Ref.  8  flow-field 

math  model,  deg  or  radians 

7  Simulated  aircraft  dive  angle,  positive  for  decreasing  altitude,  deg 

6  Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 

Xp-Yp  plane,  positive  when  the  store  nose  is  raised  as  seen  by  the 
pilot,  deg  (or  radians) 
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dp  Rotation  angle  of  the  aircraft  (flight-axis  system)  during  a 

pullup/pushover  maneuver,  deg  (see  Fig.  12) 

Ad  Angle  between  the  store  longitudinal  axis  and  its  projection  in  the 

Xp-Yp  plane,  positive  when  the  store  nose  is  raised  as  seen  by  the 
pilot,  deg 

X  Model  scale  factor 

Free-stream  density,  slugs/ft3 

<t>  Angle  between  the  store  lateral  (YB)  axis  and  the  intersection  of  the 

Yb-Zb  and  Xp-Yp  planes,  positive  for  clockwise  rotation  when 
looking  upstream,  deg 

i p  Angle  between  the  projection  of  the  store  longitudinal  axis  in  the  Xp- 

Yp  plane  and  the  XF  axis,  positive  when  the  store  nose  is  to  the  right 
as  seen  by  the  pilot,  deg 

Angular  velocity  of  the  body-axis  system,  rad/sec 

co  Angular  velocity  of  the  body,  rad/sec 

d(*)/dt  Total  derivative  of  a  vector  parameter  with  respect  to  time  such  that 

d(*)/dt  =  <5(VSt  +  d  x  C) 

-> 

<5(  )/<5t  Partial  derivative  of  a  vector  parameter  with  respect  to  time  when 

the  unit  vectors  have  fixed  directions,  i.e.,  5(  )/<5t  =  i  ( )x  +  j(  )Y 

+  kf)z 

82(*)/8t2  Second  partial  derivative  of  a  vector  parameter  with  respect  to  time 

when  the  unit  vectors  have  fixed  directions,  i.e.,  82C)/8t2  =  1  (  ”)x  + 
j  (  )Y  +  k(  )z 

(')  A  single  dot  denotes  the  first  derivative  of  a  scalar  parameter  with 

respect  to  time 
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(") 

A  double  dot  denotes  a  second  derivative  of  a  scalar  parameter  with 
respect  to  time 

0 

The  arrow  denotes  a  vector  quantity 

O' 

Primed  quantities  refer  to  model-scale  dimensions  or  properties 
defined  by  the  scaling  laws  or  to  wind  tunnel  operating  parameters 
for  the  dynamic  drop  technique 

INERTIAL-AXIS  SYSTEM  DEFINITIONS 
Coordinate  Directions 

X, 

Y, 

Z, 

Origin 

The  inertial-axis  system  origin  is  coincident  with  the  store  eg  at  release  and  translates 
along  the  initial  aircraft  flight  path  direction  at  the  free-stream  velocity.  The  coordinate  axes 
do  not  rotate  with  respect  to  the  initial  aircraft  flight  path  direction. 

FLIGHT-AXIS  SYSTEM  DEFINITIONS 

Coordinate  Directions 


Parallel  to  the  aircraft  flight  path  direction  at  store  release,  positive 
forward  as  seen  by  the  pilot 

Perpendicular  to  the  X[  and  Zi  directions,  positive  to  the  right  as  seen 
by  the  pilot 

Parallel  to  the  aircraft  plane  of  symmetry  and  perpendicular  to  the 
aircraft  flight  path  direction  at  store  release,  positive  downward  as 
seen  by  the  pilot 


XF  Parallel  to  the  current  aircraft  flight  path  direction,  positive  forward 

as  seen  by  the  pilot 

Yf  Perpendicular  to  the  XF  and  ZF  directions,  positive  to  the  right  as 

seen  by  the  pilot 
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ZF  Parallel  to  the  aircraft  plane  of  symmetry  and  perpendicular  to  the 

current  aircraft  flight  path  direction,  positive  downward  as  seen  by 
the  pilot 

Origin 

The  flight-axis  system  origin  is  coincident  with  the  store  eg  at  release.  The  origin  is  fixed 
with  respect  to  the  aircraft  and  thus  translates  along  the  current  aircraft  flight  path  at  the 
free-stream  velocity.  The  coordinate  axes  rotate  to  maintain  alignment  of  the  XF  axis  with 
the  current  aircraft  flight  path  direction. 

PYLON-AXIS  SYSTEM  DEFINITIONS 
Coordinate  Directions 


Xp  Parallel  to  the  store  longitudinal  axis  at  release  and  at  constant 

angular  orientation  with  respect  to  the  current  aircraft  flight  path 
direction,  positive  forward  as  seen  by  the  pilot 

Yp  Perpendicular  to  the  Xp  direction  and  parallel  to  the  Xf-Yf  plane, 

positive  to  the  right  as  seen  by  the  pilot 

Zp  Perpendicular  to  the  Xp  and  Yp  directions,  positive  downward  as 

seen  by  the  pilot 

Origin 

The  pylon-axis  system  origin  is  coincident  with  the  flight-axis  system  origin  and  the  store 
eg  at  release.  It  is  fixed  with  respect  to  the  aircraft  and  thus  translates  along  the  current 
aircraft  flight  path  at  the  free-stream  velocity.  The  coordinate  axes  rotate  to  maintain 
constant  angular  orientation  with  respect  to  the  current  aircraft  flight  path  direction. 

STORE  BOD  Y-AXIS  SYSTEM  DEFINITIONS 

Coordinate  Directions 


XB 


Parallel  to  the  store  longitudinal  axis,  positive  direction  is  upstream 
at  store  release 
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Yb  Perpendicular  to  XB  and  ZB  directions,  positive  to  the  right  looking 

upstream  when  the  store  and  aircraft  are  at  zero  yaw  and  roll  angles 

ZB  Perpendicular  to  the  XB  direction  and  parallel  to  the  aircraft  plane  of 

symmetry  when  the  store  and  aircraft  are  at  zero  yaw  and  roll  angles, 
positive  downward  as  seen  by  the  pilot  when  the  store  is  at  zero  pitch 
and  roll  angles 

Origin 

The  store  body-axis  system  origin  is  coincident  with  the  store  eg  at  all  times.  The  XB,  YB, 
and  ZB  coordinate  axes  rotate  with  the  store  in  pitch,  yaw,  and  roll  so  that  mass  moments  of 
inertia  about  the  three  axes  are  not  time-varying  quantities. 
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